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q 1 ABSTRACT 

<d : 

£^ ■ B02 18+357 is one of the most promising systems to determine the Hubble constant from 

time-delays in gravitational lenses. Consisting of two bright images, which are well resolved 
in VLBI observations, plus one of the most richly structured Einstein rings, it potentially 
provides better constraints for the mass model than most other systems. The main problem 
I \ left until now was the very poorly determined position of the lensing galaxy. 

After presenting detailed results from classical lens modelling, we apply our improved 
' version of the LensClean algorithm which for the first time utilizes the beautiful Einstein 

\Q \ ring for lens modelling purposes. The primary result using isothermal lens models is a now 

CN| i very well defined lens position of (255 ±6, 119 ±4) mas relative to the A image, which allows 

C^l ' the first reliable measurement of the Hubble constant from the time-delay of this system. The 

; resultof H = (78 ± 6) kms" 1 Mpc" 1 (2 cr) is very high compared with other lenses. It is, 

however, compatible with local estimates from the HST key project and with WMAP results, 
but less prone to systematic errors. We furthermore discuss possible changes of these results 
^ c*j for different radial mass profiles and find that the final values cannot be very different from the 

isothermal expectations. The power-law exponent of the potential is constrained by VLBI data 
.1 | of the compact images and the inner jet to be f3 — 1.04 ± 0.02, which confirms that the mass 

distribution is approximately isothermal (corresponding to (3 — 1), but slightly shallower. The 
' effect on Hq is reduced from the expected 4 per cent decrease by an estimate shift of the best 

, galaxy position of ca. 4 mas to at most 2 per cent. 

!• ■ Maps of the unlensed source plane produced from the best LensClean brightness model 

, ' show a typical jet structure and allow us to identify the parts which are distorted by the lens 

, to produce the radio ring. We also present a composite map which for the first time shows 

5_i ■ the rich structure of B02 18+357 on scales ranging from milli-arcseconds to arcseconds, both 

in the image plane and in the reconstructed source plane. Finally we use a comparison of 
observations at different frequencies to investigate the question of possible weakening of one 
of the images by propagation effects and/or source shifts with frequency. The data clearly 
favour the model of significant 'extinction' without noticeable source position shifts. 

The technical details of our variant of the LensClean method are presented in the ac- 
companying Paper I. 

Key words: quasars: individual: JVAS B02 18+357 - gravitational lensing - distance scale - 
techniques: interferometric 



1 INTRODUCTION 

The lensed B0218+357 jPatnaik et aill 19931) is one of the rapidly 
growing but still small number of gravitationa l lens systems with 
an accurately known time-delay l Bi ggs et alJll999t ICohen et alJ 
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2000). It is therefore a candidate to apply Refsd al's method to de- 
termine the Hubble constant Hq iRefsdal 1964). Besides the time- 
delay and an at least partial knowledge of the other cosmological 
parameters, only good mass models of the lens are needed to ac- 
complish this task. No other complicated and potentially incom- 
pletely known astrophysics enters the calculation or contributes to 
the errors. Since the cosmological parameters are believed to be 
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known with sufficient accuracy now, e.g. from the WMAP project 
iSper eel et al.Jl20 03). the only significant source of errors lies in 
the mass models themselves. The difficulties in constraining the 
mass distribution of lenses should not be underestimated, but the 
lens method is still much simpler to apply than classical distance- 
ladder methods, where errors difficult to estimate can enter at each 
of the numerous steps. 

The classical methods still fight with the problem of results 
which are i ncompatible within the forma l error bars. The HST 
key project (Mould et al.. 2000; Freedma n et al.Jl2001f) e.g. ob- 
tains a value of Ho = (71 ± 6)kms~ 1 Mpc~ 1 (la) which is 
often used as a reference. It should not be forgotten, however, that 
other groups obtained results which are significantly differe nt, even 
though the general distance ladder used i s very similar. [Sa ndage 
( 1999) determines a value of 53 ± 7 (1 a), Paro di et alJ l2000h ob- 
tain 59±6 (90 %). The differences in the analyses are various, rang- 
ing from the correction of selection effects and the Cepheid period- 
luminosity-relation to instrumental effects. The fact that a number 
of world-expert groups spent such an amount of effort but still do 
not agree on the final result proves how difficult the problems of 
distance ladder methods really are. See the references above for a 
discussi on of possible reasons for the discrepancies. In addition we 
refer to lTammann et alj 1200 2) for an extensive discussion of this 
subject and an additional result of 60 ± 5. 

After finishing a first version of this paper, the analysis of 
the results of the first year of WMAP observations have been pub- 
lished. Although Ho is not the primary goal of CMB observations, 
determinations of this parameter are possible with additional as- 
sumptions. Assuming an exactly spatially flat universe (Qtot = 1) 
and a cosmological constant for the 'dark energy' (equation of state 
parameter w = — 1), ISpergeletal]<2003l) obtain a value of 72 ± 5. 
For different w, the result for Ho can be significantly smaller. The 
assumptions of w = 1 and fi t ot = 1 can only be dropped if ad- 
ditional information from other astronomical fields are included. 
With the local information from the 2dF Galaxy Redshift Survey 
IPercival et al.Jl200lb IVerde et al.]E o02). a constraint on w is ob- 
tained. WMAP alone can constrain fi t ot only weakly, leading to 
wide allowed ranges for the density of dark energy or matter. For 
a determination of all free parameters, more astro nomical data has 
to be included as discussed in detail bv lSpergel et alJ 1200 3l) . The 
constraints for all cosmological parameters obtained from the com- 
bination of many data sets are very impressive, but it should be 
kept in mind that the situation now becomes similar to the distance 
ladder methods where model assumptions in many astrophysical 
fields are necessary, making realistic error estimates of the final re- 
sult extremely difficult. Therefore the more parameters that go into 
CMB models have accurate independent determination, the better. 
This allows the CMB analyses to concentrate on other parameters 
which cannot be determined in any other way. 

Taking these difficulties into account, we believe that the lens 
method is superior to most other methods not only in terms of the 
possible accuracy of Ho but also, and not less important, in terms 
of estimating the possible errors reliably. We emphasize that the 
effort needed is small when compared with large projects like the 
HST key project or WMAP plus 2dF survey etc. 

The only serious problem when using the lens effect to de- 
termine the Hubble constant lies in the mass models for the lenses 
which directly affect the resulting Ho- To be able to constrain a 
mass model accurately, two conditions should be met. First, the 
lens should be simple by which we mean the galaxy should be nor- 
mal without nearby companions and without a surrounding cluster. 
Only in these cases can the mass distribution be described with 



confidence by a small number of parameters. Second, as many ob- 
servational constraints as possible should be provided by the sys- 
tem. Advantageous is a high number of highly structured images. 
Extended lensed sources and Einstein rings can provide constraints 
for a wide range of positions and are therefore preferred to point- 
like multiple images. In the case of more than two images, multiple 
time-delays can also be used as constraints for the mass model. In 
addition to these pure lensing constraints, additional information 
from stellar dynamics or surface photometry of the lens can be in- 
cluded in the modelling. 

Radio loud lensed sources have the potential to provide a 
whole set of independent observational constraints additional to 
those provided by optical observations. Flux ratios of the images 
are less affected by microlensing and extinction at radio wave- 
lengths than optical ones. Radio observations also have the advan- 
tage that they can reach much higher levels of resolution than in the 
optical. Much smaller substructures in the images can therefore be 
detected and used for the lens models. 

The lens B0 2 18+357, discovered in the JVAS survey 
IPatnaik et alll993l) . meets most of the criteria for a 'Golden Lens' 
which allows an accurate determination of Hq. The lens is an iso- 
lated single galaxy, and no field galaxies nearby contribute to the 
lensing potential significantly iLehar et al. 12 000). The lens shows 
only two images of the core of the source, but substructure of the 
images can be used as further constraint. The main bonus is the 
presence of an Einstein rin g which sh ows a lot of substructure in 
detailed radio maps teiggs et all200lh . This paper deals with uti- 
lizing the information from the ring to constrain the mass models 
and to determine Hq. 

B02 18+357 has only one disadvantage: It is the system with 
the smallest image separation known, which makes useful direct 
measurements of the lensing galaxy's position relative to the lensed 
images extremely difficult. Since this parameter is of fundamental 
importance for the determination of H p, no reliable results have 
been possible so far. iBiggs etafl 1 19991) present a galaxy position 
which may be close to the correct value, but their error bars ar e 
highly underestimated as was already stated bv lLeharetaTl f2000). 

In this paper we will first discuss 'classical' lens models which 
only use parametrized data of the compact images to constrain the 
models. We will see that with this approach the constraints are not 
sufficient to determine the galaxy position with any accuracy. To 
be able to determine the Hubble constant, other information has to 
be included which is naturally provided by the highly structured 
radio ring. The best metho d available to u tilize this information is 
LensClean iKochanek & Naravanl 19921) . Because the method in 
its original form has serious shortcomings which prevent its use 
in a system like B02 18+357, where the dynamic range between 
the compact images and the ring is very high, we had to improve 
the algorithm con siderably. The LensClean method itself is dis- 
cussed in Paper I l Wucknitz 2003) while the results are presented 
here (Paper II). 

We will use the improved version of LensClean to deter- 
mine the galaxy position for isothermal models with an accuracy 
that is sufficient to achieve a result for the Hubble constant which 
is competitive with other methods but avoids their possible system- 
atic errors. Possible deviations from isothermal mass distributions, 
which are the main source of systematic error when combining the 
results from many lenses, will be discussed by presenting a pre- 
liminary analysis of VLBI observations of the substructure of the 
images. These will allow us to constrain the radial mass profile for 
power-law models. We will learn that deviations from isothermal 
mass distributions do not play a major role in B02 18+357, because 
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the constraints are already quite strong and in addition the effect on 
Ho is in this special case, where the galaxy position is determined 
indirectly with LensClean, much weaker than in most other sys- 
tems. We nevertheless discuss future work which will improve the 
results even more. 

As a secondary result we will use our best lens models and 
newly developed methods to produce maps of the source plane as 
it would look like without the action of the lens. We will show the 
source as well as its appearance in the lensed image plane on scales 
from milli-arcseconds to arcseconds. 

Finally we use LensClean to shed some light on the ques- 
tion of scattering induced effective 'extinction' in the A component 
of B02 18+357 by comparing results for data sets of different fre- 
quencies. 

This Pape r II as we l l as Pa per I are condensed versions of 
major parts of IWucknitd |2002b). For more details the reader is 
referred to that work. 



2 MASS MODELS 

In the absence of an accurate position of the lens relative to the 
position of the lensed images there is no need to use very compli- 
cated and general lens models. The effect of the lens position shows 
its full impact already with simple standard models. We therefore 
restrict most of our work to isothermal elliptical mass models. To 
introduce the notation, we write the lens equation as 



•«(*) 



(1) 



The true source position z s is the image position z shifted by the 
apparent deflection angle a.. In static single-plane lenses, the latter 
can be written as gradient of a potential tp: 



a(z) = Vip(z) 



(2) 



We use this potential to parametrize the mass models. General 
power-law models have 

i>{z) = r f> F(0) . (3) 

The radial power-law index is (3 — 1 in the isothermal case. 
The arbitrary azimuthal function F(8) describes the asymme- 
try. Elliptical mass distributions can be described by the f or- 
malism presented in lKorm ann. Schneider & Bartelmannl J 1994) or 



Kassiola & Kovner 1 1993). We prefer the approximation of ellipti- 
cal potentials l Kassiola & Kovner 1993), which have practical ad- 
vantages over the true elliptical mass distributions, because the in- 
version of the lens equation is possible analytically for elliptical 
potentials. For sufficiently small ellipticities, both approaches are 
equivalent. A singular elliptical power-law potential ('SIEP' in the 
isothermal case) is given by 



,2-/3 



to be very small. lLehar et al .12000) estimate them to be of the order 
0.01 to 0.02, which is negligible in our context. 

The time-delay Atij between images i and j is generally re- 
lated to the potential and the apparent deflection angle in the fol- 
lowing way: 



1 + Zd dd d s 



Ti = -|a(*. 



|2 



ATi. 



)r-#*o 



(6) 
(7) 



Here Zd is the redshift of the lens and dd, d s and dds are the nor- 
malized angular size distances to the lens, source and from lens to 
source, respectively. 

Isothermal model s have an interestin g analytical property, first 
described by Witt. Mao & Keeton 1 2000). For fixed positions of the 
images relative to the mass centre, the time-delay for given Ho is 
uniquely determined and does not depend on any other parameters 
of the model, as long as the images are fitted exactly. In our case 
this means that, given the time-delay, Ho for isothermal models can 
be calculated directly from the lens position zq using the equation 



Ti = 



1, 



zo 



(8) 



The dependence of AT on zq is linear. 

These simple relations show the importance of measuring the 
lens position directly. Once this position is known with sufficient 
accuracy, no explicit model fits are necessary, but the Hubble con- 
stant can be derived directly from the positions of the images and 
the lens. 

The resulting Hubble constant for isotherm al mod els using 

the time-delay of At = 10 d 5 ± d 4 l lBiggs et alJll999h is shown 

in Fig. Q as a function of zq. These values are calculated for a 

low-density flat universe (Q = 0.3, A = 0.7). They would be 

smaller by 5.5 per cent for an Einstein-de S itter universe (EdS). 

We use a lens redshift of z A = 0.6847 iBrowne et all 1 19931 
■ III ^^^^^"^^^^^^ ^^"^^^ 

Carilli, Ruoen & Yannv 1993; and others) and a source redshift of 

z s — 0.944 l Cohen, Lawrence & Blandford 2002). The older value 
of z s = 0.96 tLawren c3l99rl) would lead to a 4.2 per cent smaller 
Hubble constant. Included in the plot are estimates for the lens posi- 
tion fr om different sources. The optical positions from lLehar et alJ 
( 2000) are clearly not accurate enough. The centre of the ring is rel- 
atively well defined, but it can deviate from the position of the mass 
centre by 0'.'03 or even more, depending on the ellipticity of the lens 
and the size and shape of the source. We will see later that the error 
bars of Big gs et alJll99 9) are underestimated significantly so that 
their result for the lens position can not be used. 

The alternati ve time-delay of 10 d l ± l d 5 (2 a) from 
ICohen et alJ fcOOOl) is compatible with the result from lBiggs et alJ 
1 19991) which we used for our work because of the smaller error 
bars. The result for Ho obtained from the lCohen e t al. ( 2000) time- 
delay would be 4 per cent higher. 



•r.{z) 



(4) 



with 



+ 



(5) 



For small ellipticities e and moderate f3, this model corresponds 
to an elliptical mass distribution with axial ratio q = 1 — 2e(2 + 
P) //3, or q = 1 — 6e for isothermal models. We do not include 
external perturbations, because the external shear and convergence 
introduced by field galaxies and large scale structure are expected 



3 CLASSICAL ISOTHERMAL LENS MODELS 

Our first modelling attempts for B02 18+357 followed the classical 
route of using only parameters of the two compact components as 
constr aints. This is the same approach as presented in Bi ggs et all 
(1999). We will learn that it is not possible to determine uniquely 
the lens position with this information. The results shown in this 
section were obtained using isothermal models because the (small, 
see below) deviations from isothermality are not of importance 
as long as the galaxy position is not known with high accuracy. 
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Figure 1. Hq in units of kms 1 Mpc 1 as a function of galaxy position zq (for the concordance model, Q = 0.3, A = 0.3). The dashed line connects 
components A (at the origin) and B. Show n are also different estimates fo r zg. NIC1 and NIC 2 are optical positions from lLehar et alJ (2000). At the upper 
left are the results from fBiggs et alJ jl999l as well as the centre of the ring i Patnaik et alll993l) . Our final results with LENSCLEAN are shown by the ellipses 
(LC) to the upper left of the 'Biggs' error cross. These ellipses mark the formal 2 a uncertainties for the two parameters xq and yo with natural (smaller) and 
uniform (larger ellipse) weighting from Fig.|4] The 2 cr region for the one parameter Hq would be slightly smaller. 



Once this parameter is determined with LENSCLEAN, constraints 
on these deviations as well as the impact on Ho will be discussed 
in Sec.Hl 



3.1 Observational data 

All coordinates in this paper are measured with respect to the A 
image. The x coordinate is measured eastwards (increasing right 
ascension), y is measured northwards. 

Table Q shows a compilation of known relative positions of 
the A and B images as well as of the subcomponents 1 and 2 which 
were revealed by VLBA observations IPatnaik. Porcas & Browne) 
1995). There are significant differences in the given positions, 
but they are too small to affect the results seriously. For the 
B— A separation, we used the 15 GHz VLBA positions for the 
core components Al and Bl with their formal accuracy. The rel- 
ative positions of the ™£>components offer an additional possibil- 
ity to constrain the lens mod els. For the modelling we assumed 
an accuracy of 0. 1 mas as in iBiggs et alJ ll999l) and made inde- 
pendent calculations for the data from Ipatnaik et alJ J 19951) and 
iKemball. Patnaik & Porcas (2001). The difference between the two 
sets of positions is of the same order a s the assumed err or bar. 

The new VLBI data presented in Big gs et ail 120031) show, in 
addition to the already known two close subcomponents, the in- 
ner jet in both images with subcomponents at separations of up to 
ca. 10 mas which could potentially provide additional constraints. 
These data are not used for our classical modelling because the ac- 
curacy of parametrized model fits to the additional components are 
influenced significantly by the underlying smooth emission of the 
jet. The true errors of such fits are expected to be much higher than 
formal estimates so that we decided to exploit this new data set only 
later with LENSCLEAN. 

The magnification ratio of the two compact images is more 
difficult to obtain. Several effects can influence the results. One 
problem is scattering in the lensing galaxy which seems to decrease 
the flux density of A. This effect is stronger for lower frequen- 
cies and (with regard to the fluxes) thought to be almost absent at 



Table 1. Relative positions of components and subcomponents. 



data set 




x [mas] 


V [mas] 


(sub)comp 


5 GHz 


EVN 


a 


308.5 


130.3 


B-A 


8.4 GHz 


VLBI 


b 


309.00 


127.30 


Bl-Al 








309.32 


126.37 


B2-A2 


15 GHz 


VLBA 


c 


309.2 


127.4 


Bl-Al 








309.6 


126.6 


B2-A2 


15 GHz 


VLA 


d 


310.56 


127.11 


B-A 


optical 


HST 


e 


307 


126 


B-A 


8.4 GHz 


VLBI 


1 


1.18 


0.87 


A2-A1 








1.50 


-0.06 


B2-B1 


15 GHz 


VLBA 9 


1.072 


0.868 


A2-A1 








1.470 


0.000 


B2-B1 



a fro mlPatnaik et al]<1993l), given to 0.6 mas accuracy 

from Kemball et al 1 1200 ill , accuracy 0.09 mas in each coordinate 
c from Patnaik et al. 1995), given to 0.6 mas accuracy 
^ own fit (uniform weighting), same data as for LENSCLEAN 
e from lLehar etaflfaOOd) 

/ from Kemball etaflfcOOll) . accuracy ca. 0.05 mas in each coordinate 
9 from Patnaik et al. 1 1995), accuracy not given 

15 GHz. Secondly, flux densities of compact images embedded in 
a smooth surface brightness background cannot be determined un- 
ambiguously. We always have to expect an uncertainty of about the 
background surface brightness integrated over the beam. This error 
should be smaller than 10 per cent for the lower resolution data sets 
at lower frequencies and much less than that in other cases. Finally 
we observe the source at different epochs in the two images because 
of the time-delay. From measured light curves and the known time- 
delay I Big gs et alll999l) . the typical error of this effect is estimated 
to be less than 5 per cent. 

Table [2] shows flux density ratios for different frequencies, 
epochs and resolutions. There seems to be a systematic increase 
in the ratio with increasing frequency (discounting optical), which 
hints at scattering. Also the spread in values for 8.4 GHz and above 
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Table 2. Flux density ratio A/B for different frequencies and instruments. 



data set 




epoch 




1.7 GHz 


VLBI 


a 


1 0Q^ Inn 1 Q 


2 62 


V I rl / 




b 


1 OQ 1 Ann 

Lyy I rtUg ZO 


O OVA 

z.y id 


5 GHz 


MFRT T1\T 


b 


1 nAO t 1 1 

1992 Jan 13 


3.23 


5 GHz 




b 


1992 Mar 27 


3.35 


5 GHz 


EVN 


b 


1990 Nov 19 


3.185 


8.4 GHz 


VLA 


b 


1991 Aug 1 


3.247 


o . 4 un/ 


V L^Pl 


c 


1996/1997 


3.57 ±0.01 


Q /I /" LI ., 
0.4 <jrlZ 


\n a 
Vl^A 


d 


1996/1997 


3.2 ±0.35 


8.4 GHz 


VLBI 


e 


ivy j iviay y 


Q 1 Q 4-0 1 7 






/ 






15 GHz 


VLA 


b 


1991 Aug 1 


3.690 


15 GHz 


VLA 


9 


1992 Nov 18 


3.79 


15 GHz 


VLA 


h 


1996/1997 


3.73 ±0.01 


15 GHz 


VLA 


d 


1996/1997 


4.3 ±0.65 


15 GHz 


VLBA 


i 


1994 Oct 3 


3.623±0.065 


22 GHz 


VLA 


b 


1991 Aug 1 


3.636 


5550 A 


HST 


3 




0.14 



a from Patnaik & Porcas N 1999ft 
b from lPatnaik et alJll993l) 

c from lBiggs et alj ll99S3t) - simultaneously fitted with time-delay, used 
for their models 

from ICohen et alJ (2000), varying part simultaneously fitted with 
time-delay 

e subcomponents Al/Bl iKemball et all200ll) 

J subcomponents A2/B2 i Kem ball et all200ll) 

9 own fit (uniform weighting), same data as for LENSCLEAN 

^ from lBiggs et alJ ll999h . simultaneously fitted with time-delay 

* from lPatnaik et alj|l995t) 

5 from i'Koch anek et all2002llLehar et all2000l) 

is relatively small. For both these reasons, and because there are 
values with the effects of delay removed, we chose to use a value of 
3.75 which is compatible with most of the measurements at higher 
frequencies. Please note that the results depend only very weakly 
on small changes of this ratio. For the LENS CLEAN work presented 
below, no explicit measurement of the flux density ratio is required. 
Resolution effects and contamination by the ring are then taken into 
account automatically. 

The subcomponents are only marginally resolved and even 
the formal uncertainties of fitted Gaussians are too high to serve 
as additional constraints. Furthermore, the apparent sizes and 
shapes might be influenced by scatter broadening in the inter- 
stellar med ium of the lensing galaxy, which is known to be gas 
rich (Wiklind & Combeslll995l) . We nevertheless tried to include 
the shapes of the subcomponents into the modelling, without any 
gain in accur a cy of Hp. The actual obse r vation al data as given by 
Patna ik et all < 19951) and lKemball et all <200ll) are shown in Ta- 
ble^] The fit of the best source shapes for each lens model is per- 
formed analytically with a Cartesian parametrization of the ellipses. 
Mathematical details can be found in Appendix lA2l 

3.2 Results 

Before proceeding with numerical models, we want to compare 
the numbers of parameters and constraints to learn what can be 
achieved in the absence of degeneracies. The SIEP models intro- 
duced above have a total of five free parameters (mass scale an, 



Table 3. Shapes of the subcomponents 







major axis [mas] 


minor axis [mas] 


p.a. [dog] 


Al 


a 


0.58 ±0.05 


0.28 ±0.05 


-37 ± 5 




b 


1.90 ± 0.30 


1.12 ± 0.30 


157 ± 6 


A2 


a 


1.02 ±0.05 


0.54 ±0.05 


-47 ± 5 




b 


^ 2.0 


> 2.0 




Bl 


a 


0.36 ±0.05 


0.16 ±0.05 


-65 ±15 




b 


0.72 ±0.30 


0.35 ±0.04 


73 ±12 


B2 


a 


0.61 ±0.05 


0.23 ±0.05 


83 ±10 




b 


^ 0.73 


> 0.73 





a 15 GHz IPatnaik et alll995h 
b 8.4 GHz iKemball et all20()Ih 

position xo, yo and ellipticity e x , e y ), not counting the source po- 
sition. The positions of A and B give us two constraints (or four 
constraints minus two parameters for the true source position), the 
flux density ratio another one. The relative positions of the sub- 
components provide another two totalling in five constraints. These 
numbers may seem promising but we have to be aware of possible 
degeneracies which might prevent us from determining all param- 
eters and Ho- This is in fact the case. 

Without using the subcomponents, only three constraints are 
available. From this small number of constraints, it is impossible to 
determine the parameters of our mass model. It is in fact possible 
to find exactly one fitting model for each given galaxy position zo 
plotted in Fig.Q A wide range of (positive and formally negative) 
values of Ho is compatible with these data. 

In the next step, we include relative positions of the subcom- 
ponents from Table Q See Appendix I A II for details on how these 
data are used in the modelling. Naively, we would expect these two 
more numbers to constrain the galaxy position significantly. Fig- 
ure |2|(a)+(b) shows the total x 2 of the best fitting models for a 
range of fixed galaxy positions. For these models, we used a flux 
density ratio of 3.75, but other values measured at high frequency 
would lead to similar results. The relative parity of the two images 
was fixed to be negative as it is necessary for double image lens 
models. With our fitting algorithm, this was already sufficient to 
exclude models with more than two images in those regions where 
double lenses are at all possible. 

As the plots show, the subcomponent positions provi de only 
one m ore effective constraint, confirming the result of Leh ar et alJ 
(2000). This is easily explained. In the best models shown, the 
subcomponents are located nearly exactly radially with respect to 
the mass centre. Since singular isothermal models show no depen- 
dence of the radial stretching on image positions (there is no ra- 
dial stretching at all), we lose one constraint. The only constraint 
left acts perpendicular to the line connecting A and B. This is un- 
fortunately not the direction in which Ho changes (see Fig. 0. 
For galaxy positions leading to realistic ranges of Ho (like our 
LENSCLEAN results which will be presented later), t he minimal 
X 2 is about 1.6 for the 15 GHz VLB A data JPatnaik et alJI 19951) 
and 3. 1 for the 8.4 GHz VLBI data from lKemball et alJfcOOll) 'The 
formal number of degrees of freedom is v — in this case, but the 
model degeneracies lead to an effective value of v = 1. The fact 
that the data can be fitted quite well is interesting, because it shows 

1 If the radial mass index (3 is allowed to vary freely, the lowest residuals 
come very close to zero. 
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(a) without shapes (15 GHz) (b) without shapes (8.4 GHz) 
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Figure 2. Best x 2 for isothermal models as a function of galaxy position. Subcomponent positions and shapes for the left plots are from 1 5 GHz VLBA data 
(Patnaik et al. 1995), for the right plots from 8.4 GHz VLBI data iKemball et al. 2001). Constraints for the models are the relative position of A and B, the 
flux density ratio and the positions of the VLBA subcomponents. The upper plots are without, the lower ones with the subcomponent shapes included as 
constraints. The dashed curves separate allowed regions with two images in the central parts of the plots from regions where more than two images would be 
seen. The small ellipse marked LC is the 2 a confidence region of our final LENSCLEAN result for natural weighting, see Fig.l4ltbV 



that the data are at least approximately compatible with isothermal 
models. 

We already noted that the shapes of the subcomponents are not 
a very promising candidate for further constraints. Figure[2](c)+(d) 
confirms this view. The % 2 contribution from the shapes alone is not 
sensitive for shifts in the radial direction and the total % 2 function 
does not change qualitatively. The relatively high x 2 is a hint that 
the measured shapes are not exactly the correct ones, just as one 
might expect if the shapes are affected by scattering. The ineffec- 
tiveness of shape constraints does not change when non-isothermal 
models are used. This is a result of the fact that most information 
about the relative magnification matrix is already provided by the 
magnification ratio and the radial stretching of the subcomponent 
separations. 

We conclude that the available data on the two images are 
compatible with lens positions leading to an extremely wide range 
of values of Ho and we can therefore not confirm the relatively 
small error bars given by Biggs et al. 1 1999), even if (in contrast to 
Lehar et al. 200G) the subcomponent shapes are included as con- 
straints. Without any further information about the position of the 



lensing galaxy, 'classical' lens modelling using only the compact 
components does not provide any useful estimates for Ho . 

4 LENSCLEAN 
4.1 The method 

One promising way to get more information about the mass 
distribution of the lens is to make use of the extende d emis- 
sion which shows as a ring in radio maps (e.g. iBiggs etal] 
l200lt> . Several approaches for models usi ng extended emission 
have been discussed in the literature (e.g. iKochanek et al J | 198^: 
|^chane^^^^g^a^^2llSii^£ffiiIS^h ane k: & Hewita fl996l : 
Wallington. Kochanek & Naravan 1996). They all work by, for a 
given lens model, constructing a map of the source which mini- 
mizes the deviations from observations when mapped back to the 

2 HSTI ACS observations performed to measure the galaxy position directly 
are currently being analysed. They will provide a value for this parameter 
which might be formally less accurate than our LENSCLEAN results but 
which is completely independent of the lens model. (PI: N. Jackson) 



© 0000 RAS, MNRAS 000, 000-000 



Models for the lens and source of B02 18+357 1 



lens plane. The minimal residuals themselves are then used in an 
outer loop to find the best lens model. 

In Paper I we describe the LensClean me thod which was 
first pr oposed bv lKochanek & Naravanl 1 19921) and Ellithorp e et alJ 
( 1996) together with a number of improvements and the details of 
our own implementation. The idea of LensClean as an extension 
of standard CLEAN is very simple. CLEAN builds up a radio map by 
successively collecting point sources at the peaks of the so-called 
dirty map, subtracting them from the observed data and in the end 
convolving the collection of components with a Gaussian in order 
to produce the CLEAN map. LENSCLEAN works similarly, but for 
each position in the source plane, emission has to be subtracted at 
the positions of all corresponding image positions simultaneously 
with the magnifications given by the (for now fixed) lens model. 
In this way an emission model is built which is exactly compatible 
with a given lens model and which comprises the best fit to the data 
under this assumption. An outer loop to determine the lens model 
parameters themselves is then built around the LensClean core. 
In this loop the remaining residuals of LensClean are minimized 
by varying the lens model parameters. The result is a combined 
best-fitting model for the lens mass distribution and the brightness 
distribution of the source. 

For our work we used a data set taken at the VLA in 1992 at 
15 GHz with a bandwidth of 50 MHz and a total on-source time of 
about 6 hours. The beam size is 129 x 146 mas (p.a. — 73 °) with 
natural and 86 x 88 mas 2 (p.a. —60 °) with uniform weighting. The 
thermal rms noise per beam is 55 |jJy for natural and 140 pJy for 
uniform weighting. The same data were used in the development of 
LensClean described in Paper I. More details can be found there. 

The standard pixel and field size used for the final runs is a 
512 x 512 pixels of 5 x 5mas 2 each. A loop gain of g — 0.38 
is used and 5000 iterations of our improved unbiased LensClean 
variant are performed. Smaller pixels and a smaller loop gain would 
make the residual function smoother, which helps in finding the 
minimum, but would increase the computation times on the other 
hand. Systematic deviations caused by too large pixels or too high 
g are not significant for the parameter ranges used. For the first 
ungridded step which removes most of the bright compact images, 
a loop gain of 7 = 0.95 is used. As magnification limit for the 
images we chose a value of 100. 

As explained in more detail in Paper I, we scanned the range of 
possible lens positions zq and fitted the remaining lens model pa- 
rameters for each position in order to stabilize the fits and to be able 
to determine confidence regions for the position from the final maps 
of residuals as a function of Zq. To reduce the remaining numerical 
noise, smooth polynomial functions are fitted to the residual func- 
tion R (zo) to determine the minimum and confidence regions. 

Possible changes of the flux ratio A/B as a result of the time- 
delay in combination with intrinsic variability are taken into ac- 
count with the methods described in Paper I. From the residu- 
als of the best models we learned that any ratio changes have to 
be less than 5 per cent. This is in good agreement with monitor- 
ing data covering the epoch of our observations which showed 
an almost c onstant fl ux ratio of 3. 79 ± 0.07 during that time 
ICorbett. Browne & Wilkinsonlll996l) . The measured flux ratio is 
therefore a good estimate of the true magnification ratio, almost 
unaffected by the combination of variability and time-delay. Taken 
together, these results also confirm that, even if scattering or other 
propagation effects in the interstellar medium (ISM) of the lensing 
galaxy may change the flux ratio at lower frequencies, the effect 
is not relevant at 15 GHz and can thus be neglected here. We will 
later estimate the effect at 5 GHz by comparing with our results 



for 15 GHz. The shifts of the best zq caused by flux ratio changes 
below 5 per cent are not significant and will be neglected. 

Because of the numerical difficulties of the LensClean al- 
gorithm in combination with lens models which do not allow an 
analytic inversion of the lens equation, most of our calculations are 
done for isothermal (SIEP) models without external shear. As ex- 
plained above and in Paper I, these relatively simple models are 
already sufficient to show the effect of the very important lens po- 
sition on the determination of Ho. We will discuss later that the 
true radial mass profile seems to be very close to isothermal and 
we will estimate the effect of this deviation as a perturbation of the 
SIEP model. 

4.2 Determining the lens position 

We started the analysis with the left and right handed circular po- 
larization separately to be able to use the comparison of the two as 
consistency check. LensClean with self-calibration as described 
in Paper I was used independently on the uniformly weighted LL 
and RR data set to produce final uv data sets and residual maps. 
The residuals as a function of the galaxy position (relative to A) 
are shown in Fig. [3] To produce the plot, the residual function was 
interpolated and smoothed from the original results which were on 
a relatively coarse grid consisting of ca. 200 points. 

The confidence regions are estimated from Eq. 1 1 8i in Paper I 
and the \ 2 distribution as described there. The best normalized 
residuals are 0.77 and 0.80 for LL and RR, while the expectation 
from Eqs. 1 1 61 and 1 1 71 in Paper I would be 1 ± 0.009. Our val- 
ues are significantly lower because of the self-calibration applied 
to the data. We allowed independent phase and magnitude correc- 
tion for each integration bin independently which, by the high num- 
ber of free parameters, naturally reduces the residuals considerably. 
CLEANing and self-calibrating the data with DIFMAP without 
taking into account the lens leads to the very similar residuals of 
0.78. We conclude that the best fits with LensClean have residu- 
als within the expectations but have to keep in mind that this is not a 
very sensitive test for the goodness of fit. The best galaxy positions 
from LL and RR are clearly compatible within 1 a. 

Some care is necessary when combining the two circular po- 
larizations to obtain total intensity I which is done in the next step. 
As a starting lens model for self-calibration we used the mean of the 
best galaxy positions for LL and RR, obtained the other parameters 
for this position separately for LL and RR and took the mean of the 
two. This provides a better estimate than taking the mean of the best 
models of LL and RR directly. The LensClean calculations were 
done as before but in the self-calibration possible differences in the 
calibration errors of LL and RR had to be taken into account. To 
do this we built the emission model for I (mean of LL and RR) and 
calibrated LL and RR separately with this emission model. After- 
wards the two were again combined to obtain total intensity I. This 
standard approach relies on the assumption that no circular polar- 
ization is present and that LL and RR should (apart from noise) be 
the same. With uniform weighting only positive components were 
allowed for self-calibration. For natural weighting, we started with 
3000 uniformly weighted iterations, switched to natural weighting 
and performed another 2000 iterations. This is superior to 5000 it- 
erations with natural weighting because the improved resolution of 
the uniformly weighted beam helps in avoiding positional errors 
of the first bright components. Because of this weight switching, 
negative components were allowed as well. They are essential to 
compensate for the difference of uniform and natural weighting. 

The results are shown in Fig.|4| The plot does not show the re- 
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Figure 3. Final residuals for the 15 GHz VLA data set for (a) LL and (b) RR polarization (uniform weighting) for isothermal models. The results are smoothed 
for the contour lines (confidence limits of 1, 2, 3, 4, 5-cr). The cross hair marks the respective residual minimum, the asterisks the one of the other polarization. 
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Figure 4. Final residuals for the 15 GHz VLA data set for Stokes I for isothermal models. The results are smoothed for the contour lines (confidence limits of 
1, 2, 3, 4, 5-cr). (a) uniform, (b) natural weighting. The cross hair marks the respective residual minima, the asterisk marks the minimum for the other weighting 
scheme. The arrow in (b) symbolizes the estimated shift for a non-isothermal model with (3 = 1.04 (see Sec. 15. 21 . The shape of the residual function would 
be very similar for such a model. 



suit from the fits directly but was produced by first fitting the mod- 
els for a relatively coarse grid of zo with LensClean, then fitting 
smooth polynomials to the resulting model parameters and calcu- 
lating LensClean residuals on a finer grid where the polynomial 
fits where used for the lens model parameters. This approach does 
not change the appearance of the plot but saves a large amount of 
CPU time because real LensClean fits have to be performed only 
for a much smaller number of galaxy positions than needed to pro- 
duce a smooth plot. 

The final error ellipses are also included in Fig. Q and Al- 
though uniform weighting was superior to natural weighting in the 
first stages of calibrating the data, in the end the naturally weighted 
results are slightly more accurate if we trust the formal statistical 
uncertainties. The accuracy depends on two aspects here. One is 
the resolution of the maps, which is better for uniform weighting, 
the other is the signal to noise ratio, which is optimal for natural 
weighting. The tradeoff between the two aspects depends on the 
lens models and on the data set and it is not a priori clear which 
weighting scheme is optimal for LensClean. 

We see that the LensClean lens position is located almost 
exactly in the centre of the valley of low residuals of classical 



lens model fitting data shown in Fig. |2| The results are therefore 
consistent with each other. Since the two methods use very differ- 
ent constraints (the VLBI substructure in classical and the ring in 
LensClean modelling), this agreement provides good evidence 
that the LensClean methods actually works correctly and pro- 
duces reliable results. 

The statistical uncertainties in the total intensity (I) result 
are only slightly smaller than in the separate LL and RR re- 
sults. Slightly unexpected, the best lens position from I (uniformly 
weighted) is not located between the best positions from LL and 
RR but has an offset of ca. 6 mas, still compatible within the 1 a 
bounds. This effect is not fully understood yet. It might be related to 
low-level circular polarization in the source which would demand a 
different self-calibration technique. More probable is the presence 
of subtle instrumental effects which could not be corrected for by 
self-calibration. 

The final residuals for the naturally weighted data set and the 
best lens model are reduced x 2 = 0.842. This is again significantly 
below unity as a result of self-calibration. The value is compatible 
with results using simulated data sets but should not be used as 
test for the goodness of the fit. For the latter purpose we show in 
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Fig. |H the uniformly weighted (to increase the resolution) image 
space residuals for the best lens model (from uniformly weighted 
data) and an alternative model that is off by 1 a (calculated for 
one parameter). For these maps we used a very high number of 
100000 iterations, many more than the 5000 used for the model 
fits. We notice that in the outer singly imaged regions, where the 
CLEANing is not constrained by the lens, most of the noise has been 
included in the models so that the residuals are much smaller than 
the rms noise of the data. For the best model shown in Fig. [5](a) 
even the inner parts, which really depend on the lens model, are 
CLEANed down to a similar level. This proves that the lens model 
can not be too far from the correct one. For the alternative model, 
on the other hand, the inner residuals are significantly higher than 
before, see Fig. |5](b). These residual maps should be compared 
with Fig.|5|(c+d) in Paper I which shows the same for a simulated 
data set using a known lens model. The residuals of our alternative 
model are of the same level as in the simulations. For the best-fitting 
model, the real residuals are even smaller which could be surprising 
at the first view but which results from the fact that the best-fitting 
lens model has by definition the smallest residuals which are thus 
always smaller than residuals for the correct lens model. We con- 
clude that the total uw-space residuals and the image space residual 
map are both consistent with a perfect fit of the lens model. 

In order to test our estimates of confidence regions in the case 
of uniform weighting, we performed Monte Carlo simulations the 
results of which are shown in Fig.|4]in Paper I, The scatter of simu- 
lated results is in good agreement with the theoretical expectations 
although these could only be obtained by using approximations for 
the expected statistical properties of the residuals. In the case of 
natural weighting, the residuals follow the standard \ 2 distribution. 



4.3 Best isothermal lens model and Hubble constant 

The best fitting lens models are given in Table @ The errors are 
mainly due to the uncertainty in zp. For a fixed zp they would be 
smaller by orders of magnitude. The real flux ratio calculated from 
the models is 3.80 ± 0.03 which is compatible with the uniformly 
weighted direct model fitting to the same data set. The agreement 
with the value of 3.73 resulting from the time-delay analysis in 
iBiggs etaljj 19991) is also very good. 

The ellipticity of the mass model is unexpectedly large. We 
have to remember that the ellipticity of an equivalent elliptical mass 
distribution is three times the ellipticity of the potential which leads 
to a best estimate of the axial ratio of a/b ~ 1.4. In contrast to 
this, the lensing galaxy appears quite round on optical images. The 
uncertainties in e are quite large, however, so the real ellipticity 
might be smaller. For ellipticities as high as observed here, it is of 
concern how far the elliptical potential still is a good a pproxima- 

/ Wif 



tion fo r elliptical mass distributions. As was shown bv Iwitt et alJ 
(2000), the magnification in self-similar isothermal models is a di- 
rect function of the surface mass density which means that the crit- 
ical curve shown in Fig.|6|is also an isodensity contour. The shape 
of this curve looks rather elliptical and does not show hints of the 
unrealistic dumbbell shapes which are observed for more extreme 
ellipticities. We conclude that the use of elliptical potentials is still 
justified in our case. 

The Hubble constant Hp can now be determined from the dis- 
tances and from the parameter AT given by the lens model. For the 
final result we only use the naturally weighted data since the un- 
certainties are smaller here. For an EdS world model, the result is 
Hq = (74.1 ± 5.5) kms~ Mpc~ . For a low-density flat standard 



Table 4. Parameters of the best fitting isothermal lens models for uniform 
and natural weighting (final result VLA 15 GHz Stokes I). The error bars 
are statistical 2 a limits for one parameter each; they originate mostly from 
the uncertainty of the lens position and would be much smaller for fixed Zo ■ 
The position angle of the major axis 9 is measured in the usual astronomical 
sense from north through east. 



parameter 



uniform 



natural 



xo 


[arcsec] 


0.2582 


0158 
0087 


0.2548 


yo 


[arcsec] 


0.1210 ±0 


0057 
0048 


0.1191 


eta 


[arcsec] 


0.1616 If, 


0039 
0124 


0.1634 


e x 




+0.0043 to 


0123 
0438 


+0.0091 


ey 
M 




-0.0707 to 
0.0709 


0273 
0637 


-0.0585 
0.0592 


e 


[deg] 


-46.7 




-49.5 


AT 


[arcsec 2 ] 


0.0393 to 


0055 
0032 


0.0380 



hO.0031 
-0.0039 



t-0.0026 
-0.0034 



W.0196 
-0.0221 



model with Q — 0.3 and A = 0.7, it increases by about 6 per cent:' 



Ho = (78.5±5.8)kms _1 Mpc _ 



(9) 



The error bars are formal 2 a limits which include the uncertainty 
of the lens model (6 per cent, contributed almost exclusively by the 
uncertainty of the lens position) and the time-delay (4 per cent). 
We stress that they do not include possible systematic errors which 
might result if the true mass distribution can not be described by our 
isothermal models (SIEP). As we show below, no large effects from 
this are expected. Using the time-delay from lCohen et alJ (2000). 
our result becomes Ho = (82±13)kms~ 1 Mpc -1 with the largest 
error contribution (15 per cent) coming f rom the time-delay. 

Standard quintessence models (e.g. lLindeJl988albl) have pa- 
rameters between the EdS and the low-density A world models and 
would lead to Hubble constants between the two extremes. For 
w = —1/3, the result would be indistinguishable from the EdS 
case (equivalent to w = 0) while for w = —2/3 the Hubble con- 
stant would be 2.1 per cent higher than for EdS. The cosmological 
constant corresponds torn = —1. 



5 NON-ISOTHERMAL LENS MODELS 

It is a well known fact the the radial mass distribution is 
the most important source of uncertainties in the Ho-At re- 
lation for lens systems where the lens position is known ac- 
curatel y (see e.g. references in the introduction of Kochanek 
l2002ah . For power-law models a simple scaling law of Hp oc 
2 — is valid in many ca ses iWambsganfi & Paczvnskilll994l 



IWitt. Mao& Schechte^ fl995lWucknitz & Refsdalll200lh while in 
quadruple lenses an even stronger dependen ce of Hp o c (2 — ) //3 
was found under certain circumstances by Wucknitz 1 2002a) for 
generalized power-law models following Eq. As a generaliza- 
tion of the SIEP models, we use elliptical power-law potentials 
which (if not too far from isothermal and not too elliptical) are a 
good approximation to elliptical power-law mass distributions. The 

3 The result from the uniformly weighted data would be larger by 3.4 per 
cent which is consistent with the naturally weighted result within 1 a. 
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Figure 5. Uniformly weighted image plane residuals after 100 000 iterations of LENSCLEAN with (a) the best lens model, (b) a lens model that is excluded by 
1 <x. The expected noise of the dirty maps is 0. 14mjy per beam. Compare these maps with their simulated equivalents in Fig.|2](c+d) of Paper I. The critical 
curve and the two bright images are marked. 



two families of models differ only in the special form of the az- 
imuthal function F(ff). 

We will see that deviations from isothermality are very small 
so that our LENSCLEAN results for isothermal models can still be 
used as a first approximation. The deviations can thus be treated as 
a small perturbation. 

5.1 Constraints from classical modelling 

The subcomponents of both images are oriented more or less ra- 
dially with respect to the expected lens position. This has the 
disadvantage that their positions cannot be used to constrain the 
galaxy position for isothermal models as shown before. For non- 
isothermal models, on the other hand, this insensitivity can be used 
to determine constraints for the radial mass profile from these data 
relatively independent of the true lens position. If this position is 
taken to be close to the LENSCLEAN result and compatible with 
the VLBI constraints, /3 can be estimated with some accuracy us- 
ing the method presen ted in Append ix I A II If we use the 8.4 GHz 
VLBI data from lKemball et alJ^OOlh . the result is B = 1.04±0.02 
(1 a limits). The error bar does not include uncertainties of the 
galaxy position but the values are r elatively i ndependent of zo any- 
way. From the 15 GHz data from lPatnaik et all ll995h. a slightly 
larger value of /3 = 1.06 ± 0.03 can be estimated. Big gs et alJ 
(2003) were able to detect the jet with many subcomponents over 
a length of 10 mas in both images with global VLBI observations 
at 8.4 GHz. A systematic relative stretching of the B jet image by 
about 10 per cent can be accounted for by a non-isothermal lens 
model with B w 1.04, in very good agreement with the other re- 
sults. We conclude that the lensing galaxy seems to be close to 
isothermal with a deviation in B of only about 4 per cent, but have 
to keep in mind that these results are only preliminary estimates 
which are not based on a self-consistent lens model but used the 
lens position from LENSCLEAN with isothermal models to esti- 
mate deviations from isothermality. The errors of this method are 
expected to be only small (see also below) but a combined fit of 



non-isothermal lens models to the VLA and VLBI data should be 
aimed for in the future. 



5.2 Preliminary results from LensClean 

As explained in Paper I, LENSCLEAN relies on a very robust 
method for the inversion of the lens equation. Using the newly de- 
veloped LENTIL method implemented into LENSCLEAN, we were 
now able to perform a few tests with non-isothermal lens mod- 
els. We repeated the calculations described before with the 15 GHz 
VLA data but used different fixed values of B and compared the 
results afterwards. We noticed that the best lens position does in- 
deed shift slightly eastwards when increasing /3, by about 1 mas for 
1 per cent change of B. See the arrow in Fig.[4]for an illustration. To 
our relief, this effect acts on Hq in the opposite direction as the scal- 
ing with 2-/3, which we expect for fixed galaxy positions, so that 
the two cancel almost completely. A rather conservative estimate is 
that the 2 — 3 scaling is reduced by at least a factor of two as a result 
of the shifting lens centre. This means that for a value of 3 = 1.04 
a change of the Hubble constant by at most 2 per cent is expected 
compared to the isothermal models. This is much smaller than the 
statistical uncertainties of the result for isothermal models and thus 
not of concern yet. When increasing 3 slightly, starting from the 
isothermal 3 = 1, the resulting shift of the best LENSCLEAN lens 
position has the tendency to decrease the optimal 3 in classical fits. 
Combined fits of both data sets should therefore be possible and 
lead to very stable results. 

In the near future we will use LENSCLEAN and LENTIL 
to obtain more accurate constraints on the radial mass profile 
from the exi sting very high-quality VLBI observations at 8.4 GHz 
(Biggs et alj|2003|) . In combination with the existing VLA data 
which are more sensitive to the galaxy position, a self-consistent 
lens model can be found and accurate error estimates will be possi- 
ble. 
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Figure 7. Source plane maps of B02 18+357. (a) Restored with the anisotropic approximated source plane CLEAN beam (see text), (b) Restored with circular 
beams (similar to the lens plane be am shown at l ower right). The contour levels are the same as in Fig.[£] Also shown are the diamond shaped caustic and the 
elliptical 'pseudo-caustic' or 'cut' iKovner 1987). 



6 SOURCE AND LENS PLANE MAPS 

In Paper I we discussed a new method to reconstruct source plane 
maps and superresolved lens plane maps from emission models de- 
termined with LensClean. We present the resulting lens plane 
maps of B02 18+357 in Fig.© One version (a) is restored with the 
nominal CLEAN beam while the other (b) is a superresolved ver- 
sion with anisotropic and varying beams calculated as explained 
in Paper I. The superresolved map is not easy to interpret because 
of the varying and anisotropic beams. The resolution is generally 
higher close to the critical curve i.e. in parts of the ring. The el- 
liptical shape of the B Image is also a result of the local elliptical 
beams. 

The reconstructed source plane map is shown in Fig. The 
left map (a) is reconstructed with the source plane beams calcu- 
lated from Eqs. 1291 and 1301 in Paper I. The radial streaks emerg- 



ing from the region near the caustic are artifacts of the varying and 
highly anisotropic beams. The elliptical shape of the core is also 
a result of the highly eccentric local beam and should not be mis- 
interpreted as a resolved elliptical core. The right map (b) was re- 
stored with varying but circular beams (the smallest circles which 
still cover the local elliptical beams completely) to make interpre- 
tations easier. Unfortunately the greatest part of the lensing magni- 
fication and improvement of resolution is lost in this map. 

We see that the jet, which by coincidence emerges from the 
inner core in a direction very similar to the major axis of the el- 
liptical core image in Fig. 0(a), bends southwards and crosses the 
caustic of the lens at the bend. It is this part of the source which 
forms the ring in the lens plane and which is used as constraint by 
LensClean. On larger scales (shown especially well by longer 
wavelength observations, e.g. the 8 GHz VLA data in Fig.[S}, the 
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jet proceeds in a southern direction only to bend eastwards again 
at a distance of about 1 arcsec. The jet then seemingly ends in a 
relatively bright blob 1.5 arcsec south-east of the core. 

The composite colour map in Fig. [8] shows B02 18+357 on dif- 
ferent scales both in the lens plane and in the source plane (re- 
constructed with LensClean). The main maps showing the outer 
parts of the jet on arcsec scales are made from a 8.4 GHz data set 
obtained at the VLA simultaneously with a global VLBI observa- 
tion iBiggs et alJ [2003. Fig. 4). We show a regular lens plane map 
and a source plane reconstruction with circular beams. The insets 
on the left are made from the 15 GHz VLA data set that was used 
for our main LensClean work. The lens and source plane maps 
show the same data as Fig. |6|(a) and0(b). We also show the crit- 
ical curve in the VLA lens plane maps and the tangential caustic 
(diamond shape) and cut (elliptical) in the VLA source plane maps. 

The insets on the right (upper: B, lower: A) are regular lens 
plane CLEAN maps of a global VLBI data set iBiggs et all 20031. 
The source plane map is a projection into the source plane of 
the A image alone. The B image is less resolved and would not 
provide much additional information. The map was restored from 
th e CLEAN compone nts with the projected CLEAN beam. Fig. 5 
in Bi ggs et all J2003) shows projections of images A and B re- 
stored with equal beams to show the similarities. True LensClean 
source plane maps of the VLBI data have not been produced yet. 

The magnifications relative to the main maps are 3 for the left 
insets and 130 for the jet insets on the right. The scales are the 
same for lens plane and source plane. The contour lines start at 
three times the rms noise and double for each new level. The sur- 
face brightness levels in the lens plane and source plane maps are 
directly comparable. The magnified sections are marked by black 
rectangles in the VLA maps. 



7 INVESTIGATION OF FREQUENCY DEPENDENT 
FLUX RATIOS 

It is a generally accepted fact that the apparent flux ratio A/B 
changes systematically with the frequency at which it is measured. 
At the high end the value seems to reach its final value of 3.8 for 
frequencies > 15 GHz. On the other end the ratio goes down to 
2.6 at 1.7 GHz (compare Table At even lower frequencies it 
is difficult to separate the compact emission from the ring which 
becomes more dominant at low frequencies. There is unambigu- 
ous evidence that the ISM of the lensing galaxy is very rich and 
that optical and radio emission are affected to a significant degree 
I Wikli nd & Combeslll995l) . lBiggs et ail 120031) estimated the scat- 
tering measure from the broadening of image A relative to B (tak- 
ing into account the magnifying effect of the lens). They obtain a 
value of SM ~ 150kpcm -20,/3 which is very high compared to 
typical lines of sight along our galaxy but comparable to lines of 
sight through the Galactic Centre. 

One possible explanation for the low flux ratios at low fre- 
quencies is therefore the action of the very rich ISM of the lensing 
galaxy in front of image A. This can cause an effective extinction of 
A either by a high amount of scattering or by other physical effects. 
The refractive index of an interstellar plasma goes proportional to 
A 2 which would explain the stronger effects for larger A. 

Another possible explanation is a very strong frequency de- 
pendence of the source structure which would effectively result in a 
position shift along the jet for lower frequencies. Because the spec- 
trum of the jet is much steeper than the spectrum of the core, such 



an effect is not implausible. When the source is shifted in this direc- 
tion, the situation becomes more symmetric and the magnification 
ratio decreases. For a sufficiently large shift (> 10 mas) the differ- 
ential magnification gradient could, in combination with this fre- 
quency dependent source position, well explain the observed flux 
ratio changes. There are, however, strong arguments against this 
scenario. If the shift is as large as required, there should be signifi- 
cant differences in the appearance of the innermost compon ents in 
VLBI maps at 8.4 GHz (see e.g.lKemball et alJ[200l|:lBiggs et all 
2003) and 15 GHz (see e.g. Patnaiketal. 1995), but these maps 
look remarkably similar. Both show two strong inner components 
with similar shapes and compatible relative positions. Their sep- 
aration of ~ 1.4 mas is far too small to explain the demanded 
shift by different spectral indices. If the shift is caused by another 
component which is not seen at 15 GHz but becomes stronger at 
low frequencies, such a component should be detectable on a scale 
of 10 to 20 mas in the e xisting very deep 8.4 GHz VLBI observa- 
tions. [Bigg^yjl] GoIi3) d° indeed detect more jet components but 
these are far too weak to account for a shift of the required magni- 
tude. A more exotic alternative scenario would be refraction caused 
by large scale systematic trends in the column density of the lens 
galaxy's ISM. 

To investigate the question of extinction and/or source shifts 
with LensClean, we compare our best lens models fitted to the 
15 GHz data which are known to be unaffected by both effects with 
a MERLIN data set at 5 GH z. This multi-frequ ency data set was 
used to produce the maps in Bi ggs et all 1200 ll) and is described 
in detail there. We used only the MERLIN part because the VLA 
resolution at 5 GHz is too low to be of any help here. 

Unfortunately, this data set has some disadvantages for 
LensClean. While the resolution is better than in the 15 GHz 
VLA data, and the flux in the ring is higher due to the lower fre- 
quency, it was not possible to obtain reliable direct constraints for 
the lens model with LensClean. One problem is the possibility 
of extinction itself. Another very serious problem is the frequency- 
dependence of the emission, which i s different for the ring and the 
compact images. IBiggs etal] 1200 ll) approximately corrected for 
this effect by first mapping the three frequencies independently. 
The CLEAN components responsible for the compact images were 
then subtracted from the data so that only the ring remains. The 
three data sets were rescaled in magnitude to obtain a consistent 
total flux density at all frequencies. Finally, the compact compo- 
nents from the central frequency were added back to the complete 
data set. This process may introduce some distortions to the low- 
est and highest frequencies, because it combines data for the ring 
from these frequencies with data for the compact components from 
the central frequency. The combined data were then successively 
mapped and self-calibrated to obtain the final data set and map. 

We used the resulting uv data set as basis for our computa- 
tions. The spectral index correction seems to work very well for 
making maps, but introduces errors in LensClean which are not 
fully understood. This shows especially in incompatible results 
from the individual three frequencies and the combination of all 
three. 

On the other hand, tests showed that the data can nevertheless 
be used to fit a subset of the lens model parameters if the other 
ones are fixed. At most three parameters can be determined ac- 
curately because they are already well defined by the two bright 
images alone and the ring is only a relatively small correction. We 
therefore used a variety of lens models compatible with the 15 GHz 
VLA data and fitted the same models to the 5 GHz data set allow- 
ing for a relative shift of the two data sets and a possible extinction 
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Figure 8. Composite image of the lens plane (top) and source plane (bottom) of B0218+357. The main central images show 8.4GHz data from the VLA (tick 
marks 1 arcsec). The insets on the left are from the 15 GHz VLA data set (tick marks 100 mas, compare Fig.|6|and0, the jets on the right are from 8.4 GHz 
global VLBI observations (tick marks 10 mas). The magnified regions are marked by black rectangles. See Sectionlolon page !12l for details. 
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Figure 9. Fits of MERLIN 5 GHz data with lens models fitting the VLA 
15 GHz data as a function of the extinction in A. top: residuals, bottom: 
relative shift of MERLIN— VLA data. 

in the A image at 5 GHz. These fits are much more robust than 
complete lens model fits from this data set alone. In particular, the 
results do not depend on which of the frequencies is used. The com- 
bination of all three also leads to the same result. 

Our strategy to determine shift and extinction works like this. 
First we took one of the best lens models from LENSCLEANing the 
15 GHz data. Then we fixed all parameters but zq and made a fit 
with the 5 GHz data including a (for the moment fixed) extinction 
in A. The difference between the two values of zo is then a mea- 
surement of the shift in the lens plane because in reality the lens 
should have the same position at both frequencies but the source 
may shift. We then varied the extinction to find the best value by 
minimizing the residuals. 

The results for one of the best VLA 15 GHz lens models are 
shown in Fig. [9] The position of the minimum changes some- 
what when other lens models still compatible with the 15 GHz data 
are used. The most extreme values for the extinction are 15 and 
35 per cent; most probable is a value between 20 and 30 per cent. 
The data sets where initially registered to obtain the same position 
for the A image. The relative shift in the fit (measured as a shift 
of the lens centre MERLIN — VLA, equivalent to a shift of the 
source by the same amount with opposite direction) in the region 
of the minimum is negligible. The fact that the shifts in x and y have 
their zeros at almost the same extinction which itself is furthermore 
compatible with the residual minimum gives strong evidence that 
the real shift is indeed small and that the brightness of A at 5 GHz 
is reduced by ~ 25 per cent by some propagation effect, possibly 
scattering. For zero extinction, the required shift would be about 
11 mas eastwards, but the residual curve seems to be incompatible 
with this scenario 4 . It seems as if extinction in A, caused by what- 
ever physical process, accounts for a least most of the frequency 
dependent flux ratio. This has an important consequence for fur- 
ther modelling work. Although the ring becomes brighter at lower 
frequencies, only high frequencies can be used in a simple way to 
fit lens models. Observations at medium frequencies may on the 

4 We do not give a formal significance or confidence limits of extinction 
or shift because the systematic problems with this data set are not fully 
understood yet. 



other hand be used to study the propagation effects in the ISM of 
the lensing galaxy in more detail now that the lens model is known 
with sufficient accuracy. 



8 DISCUSSION 

B02 18+357 is one of the most promising lens s ystems fo r deter- 
mining the Hubble constant with the method of iRefsdal 11964) . 
This method has, compared with other approaches, the advantage 
that it is a very simple one-step determination and relies on the un- 
derstanding of only very little astrophysics. In this way the system- 
atic uncertainties can not only be minimized but (equally impor- 
tant) also be estimated much better than in distance ladder methods 
which have their problems with several astrophysical processes at 
each step. The only serious uncertainty left is the mass distribution 
of the lens. In this respect B0218+357 is a close to optimal 'Golden 
Lens'. The lensing galaxy seems to be a regular isolated galaxy 
without contributions from a group or cluster nearby. This not only 
avoids the inclusions of further parameters to describe the exter- 
nal mass distribution but also allows the assumption that the mass 
distribution of the galaxy itself is that of an unperturbed smooth 
ellipsoid. B02 18+357 provides a wealth of potential constraints for 
the mass models in the form of highly accurate image positions, 
a measurable flux density ratio with probably negligible influence 
by microlensing, substructure on scales resolved by VLBI and a 
well resolved and structured Einstein ring on the scale of the image 
separation. 

Without using the ring, 'classical' lens modelling allows us 
to tightly constrain the ellipticity and mass of the galaxy as long 
as its position is known. The VLBI structure of several subcom- 
ponents along the jet in a radial direction relative to the lens al- 
lows us to determine the radial mass distribution with an accuracy 
not p ossible in most other lenses (see e.g. lKochanek & Schechteil 
2003, and references therein). This most important general degen- 
eracy of the lens method can therefore be broken in the case of 
B0218+357. The only real disadvantage is the small size (~ 0'.'3) 
of the system. B02 18+357 is indeed the lens with the smallest im- 
age separation known to date. Although the images and the lens 
are detected with optical observations, it was until now not pos- 
sible to use these observations to accurately measure the galaxy 
position directly. HST/ ACS observations obtained recently are cur- 
rently analysed to allow a first usefu l direct meas urement. 

The model lens positions from lBiggs et alJ ll999h are proba- 
bly not too far from the truth and are indeed consistent with our 
results but their error bars are seriously underestimated. We have 
shown that, without using the ring, no useful estimate of the galaxy 
position and thus Ho is possible. 

Previous efforts have not utilized the most striking feature of 
B02 18+357, the beautiful Einstein ring shown by radio observa- 
tions (e.g. IBiggs et alJl2001r) . Doing this is much more difficult 
than using parameters of compact components because the source 
itself has to be modelled as well, either implicitly or explicitly. 
The method best suited for this task is LensClean which sur- 
prisingly has been used only for very few cases before. This is 
partly a result of the high numerical demands and partly of seri- 
ou s sh ortcoming s of the original met hod as it was imple mented 
bv lKochanek & Naravanl ll992h and lEllithorpe et alj<1996h . In Pa- 
per I we discuss the development of our new improved variant of 
LensClean. One of the most important improvements is the im- 
plementation of a new unbiased selection of components. 

The problem of determining the lens position is almost de- 
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generate because the bright components, which contribute most of 
the signal, provide no information for this parameter. It is only the 
relatively weak and fuzzy ring which can be used as constraint. In 
other words the residuals have a very strong dependency on two 
directions (defined by the compact components) of the five dimen- 
sional parameter space but change only weakly in the other three 
directions as a result of the constraints provided by the ring. Finding 
minima of functions like this is a numerically difficult problem and 
relies on the very accurate knowledge of the residual function. Any 
numerical noise produces local minima which make the finding of 
the global minimum difficult. Special care was therefore needed to 
produce the best possible results with the given limited computing 
resources. 

The modified LensClean algorithm was applied to a 15 GHz 
VLA data set resulting in good constraints for all parameters of an 
isothermal lens model, including the lens position of xo ~ (255 ± 
6) mas, yo — (119±4) mas relative to the A image. This model was 
th en used to determi ne the Hubble constant from the time-delay 
of iBiggs et all Jl999h to be H = (78 ± 6) km s" 1 Mpc" 1 . The 
error bars are 2 a confidence limits which for the Hubble constant 
include the error of the lens position (6 per cent) and the time delay 
(4 per cent). The accuracy for all other lens model parameters is 
much higher as long as the position is fixed. The value for Ho is in 
agreement with the results from the HST key project I Mould et al. 
2000; Freedman et al. 2001 ) and the WMAP project l Spergel et al. 
2003) but incompatibl e with the lo we r dis tance-ladder resul ts of 
Sandal ll999h . |Parodi et alj feOOOh or lTammann et alJ |2002h. 

Recently, a series of papers was published by Kochanek 
(2002Jbl E003h . see also lKochanek & SchechteJ <2003h . in which 
it is claimed that a number of gravitational lens time-delays is 
compatible with the Hubble constant from the HST key project 
only if the mass concentration is as compact as the light distribu- 
tion (J3 < 1 in the picture of power-law models). Our result for 
B02 18+357 does not confirm this view and does not give rise to 
a 'new dark matter problem' since it is significantly higher than 
other results from lenses. Our isothermal mode ls lead to a valu e 
absolutely compatible with the one preferred in lKochanekl f2002a1. 
O ne could now in a similar way c ompare with the lower Ho results 
of lSandagd <1999h . |Parodi et alJ <200Cl) or lTammann et alJ <2002t) . 
However, as long as the local determinations do not agree with each 
other, we consider such an exercise to be of only limited value. 

Rather than using locally measured values for Ho to con- 
strain lens mass distribution, we prefer a more direct approach, 
either using the lens effect itself (see below), or by including 
additional information. The latter approach is follo wed by thy 
'Lens e s Structure and Dynamics' su rv ey LS D I KooDmans & Treu 
2002t iTreu & Koopmans! l2002biat iKoopmans & TreJ l2003a ; 
Koopmans et al. 2003; Koopmans & Treu 2003b). The general idea 
is simple: The total mass within the Einstein radius is well con- 
strained by lensing, in contrast to the radial mass profile. For the 
given mass, the stellar velocity dispersion of the lensing galaxy de- 
pends strongly on this profile. The more the mass is concentrated in 
the centre, the deeper the potential well and the higher the velocity 
dispersion has to be. With some additional assumptions, it is indeed 
possible to obtain valuable constraints which can then again be in- 
cluded in the lens models to determine Hp. Fo r the important sys- 
tem PG 1 1 15+080. lTreu & Koopmans| l2002bli obtain a power-law 
index of f3 = 0.65 which increases the otherwise very low estimate 
of Ho for this system to a value of 59 ± 10 (1 o). This steep mass 
profile is quite unusual and moves the galaxy significantly off the 
fundamental plane. Nevertheless it is not as steep as the light profile 
and the result for Ho correspondingly lies between the isothermal 



and constant M/L model results o flKochanek & Schechted 12003). 
For B 1608+656, Koopmans et al. 1 2003) are able to constrain the 
mass profile of the main lensing galaxy quite well and find that it is 
consistent with isothermal. The resulting Hubble constant is 75 ± 7 
(1 o), inconsistent with the v ery low mean value of 48 ± 3 obtained 
bv lKochanekl l 2002hLE00ll for isothermal models for a number 
of lenses, but consistent with our result for B0218+357. However, 
B 1608+656 has the major disadvantage of a second lensing galaxy 
very close to the primary. The two galaxies might even be interact- 
ing which could cause complicated deviations from usual galaxy 
mass distributions. 

These examples show that the general picture of lensing con- 
straints for Ho and m ass profiles is not as consi stent as it appeared 
not long ago (e.g. Koopmans & Fassnacht 1999). In our opinion the 
only way to resolve the discrepancies is to try and use the best avail- 
able constraints for all applicable lens systems individually. Our 
work on B02 18+357 shows this for one example, although there 
are still many open questions. 

Our alternative way of constraining mass profiles relies on 
the lens effect alone and has the advantage that additional astro- 
physics (galactic dynamics in LSD work) and corresponding addi- 
tional model assumptions are not needed, keeping the method sim- 
ple and clean. It is easy to understand that lenses showing only 
two or four images of one compact source cannot provide suffi- 
cient information to constrain the mass profile tightly, especially 
since quads usually have all images close to the Einstein ring so 
that they probe the potential at only one radius. Far better suited 
are systems with extended sources or at least with substructure in 
the compact images. B0218+357 is a good example which offers 
both. 

Using our galaxy position as an estimate, non-isothermal 
power-law models could be constrained with the VLBI substruc- 
ture. These are not the final results since no self-consistent fitting 
of all available data has been performed so far, but the weak depen- 
dency of the radial mass exponent /3 on the galaxy positions shows 
that the estimate is nevertheless quite good. We learned that the 
mass distribution in the lens is close to isothermal (which would 
be /3 = 1) but is slightly shallower {J3 ~ 1.04). Preliminary cal- 
culations with LensClean showed that the shift in the best lens 
position caused by this small deviation partly compensates for the 
expected scaling with 2-/3 for a fixed position. We therefore do 
not expect very significant changes of the result for Ho from the 
slight non-isothermality in this system. A very conservative error 
estimate disregarding this compensation effect would be a possible 
4 per cent error for Ho- 
lt has to be kept in mind, however, that all results presented so 
far depend on the assumption of a power-law for the mass profile 
(and potential and deflection angle). For more general mass models 
the constraints on /3 are still interesting, but they only measure the 
local slope 5 of the profile rather than the global mass profile. It is 
expected that power-laws can be used as a very good approximation 
for more general models for limited ranges of radii, and our model 
fits confirm this view in the case of B0218+357. The LensClean 
results for the lens position presented here are therefore also valid 
for other mass profiles, but the value for the Hubble constant may 
well have to be modified slightly. 

The mass models discussed here do not include external shear 
because it seems to be very small. The expected true external shear 



5 To be more specific the VLBI substructure provides constraints for the 
difference of the slopes of the deflection angle at the two image positions. 
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of 1-2 per cent iLehar et all200Ch would change the result of Ho 
by about the same amount at most. More difficult to estimate is the 
possible influence of differences in the ellipticity of the inner and 
outer parts of the lensing galaxy which can effectively also act as 
external shear. The possible impact of the remaining higher-order 
aspects of the mass profile degeneracies and the effects of multi- 
component mass models (bulge+disc+halo) are currently under in- 
vestigation, especially the aspect of what apparent isothermality in 
the central part of the galaxy, which should be dominated by the 
luminous matter in the bulge, means for the true mass distribution. 

In the future we will avoid the model fitting of VLBI 
components which are then used as constraints and instead use 
LensClean itself also on the existing VLBI data that show a 
wealth of structure in the images (see Fig. [S). To be able to use 
non-isothermal models with LensClean, we developed the new 
method LENTIL for the inversion of the lens equation. With these 
two methods combined, a simultaneous fit of medium resolution 
VLA data (sensitive t o the lens position ) and the existing high res- 
olution VLBI data of Bigg s et al] 120031) . which are more sensitive 
to the radial mass profile, will not only improve the results for all 
parameters but will, by using self-consistent lens models, also pro- 
vide realistic error estimates including all uncertainties. 

To be able to use the 8.4 GHz data set it will be necessary to 
include the extinction and scatter broadening in the A image in the 
models. This is possible because the 15 GHz data are only little af- 
fected by extinction so that the combination of the two allows esti- 
mates of this effect. Since the VLBI constraints for the lens models 
are mainly given by the positions of subcomponents but not by their 
size, which is relatively independent of the lens model, the different 
sizes of components in A and B can be used to obtain better con- 
strain ts on the sca ttering measure than the simple estimates from 
iBiggs etaljj2003l) . 

The application of LensClean to the VLBI data will also 
show whether any significant amount of substructure in the mass 
distribution of the lens is needed t o expla in all features of the jet. 
The first results from Big gs et all 120031) seem to be compatible 
with no substructure, but a quantitative analysis can be done in the 
future. If significant substructure effects are present, it will also be 
necessary to correct the medium resolution data for them before 
using LensClean to determine the galaxy position. Clumps of 
matter close to one of the images could change the magnification 
significantly which would then mimic a different lens model. The 
analysis of the VLBI data will show whether this may be the ex- 
planation for the (compared with some other lenses) relatively high 
result for Ho. LensClean can also be used to study substructure 
from VLBI observations of ot her lenses. It w ill then be possible to 
do an analysis similar to e.g. M etcalfl 120021) but without any as- 
sumption about the true source structure. 

To improve the medium resolution side, we recently per- 
formed observations of B0218+357 at 15 GHz using the fibre link 
between the VLA and the VLBA telescope Pie Town. This combi- 
nation provides a resolution 1 .5-2 times better than the VLA alone. 
The long track observation of the total accessible hour angle range 
of 14 hours also improves the sensitivity which in combination with 
the improved resolution can reduce both statistical and (even more 
important) systematical errors. Simulations performed before the 
observations showed that an accuracy of about 1-2 mas for the lens 
position should be achievable. This would in comparison to the ex- 
isting data be an improvement of a factor of ~ 5. The remaining 
uncertainty in Ho from the lens position alone will then shrink to 
below 2 per cent. Possible uncertainties caused by calibration errors 
will also be reduced significantly. This new data set is currently be- 



ing calibrated and will be analysed with LensClean soon, alone 
and in combination with the VLBI data. 

As a secondary result of LensClean we presented maps of 
the brightness distribution of the source as well as improved lens 
plane maps, both produced with a method newly developed to- 
gether with LensClean in Paper I. This allowed the first view 
of B0218+357's source as it would be seen without the distortion 
of the lens, but with improved resolution. The ring itself is caused 
by a bending jet which crosses the tangential caustic of the lens. 
In the future we will improve these results and also produce source 
plane polarization ma ps. The ring shows an interesting radial polar- 
ization pattern iBiggs et all2 0031 and the data will, when projected 
to the source plane, give us a very detailed and magnified view of 
the polarization structure in the jet of the source. 

Finally we used LensClean to investigate the puzzling 
changes of the flux density ratio A/B with frequency. To test the 
two theories of either scattering induced extinction caused by the 
ISM in the lensing galaxy or of an effective shift of the source with 
frequency, we compared our best LensClean models, which were 
fitted to 15 GHz data where the fluxes are expected to be unaffected 
by propagation effects, with MERLIN data at 5 GHz taking into 
account a possible relative shift and extinction in the 5 GHz data. 
Although the 5 GHz data set has some calibration problems, the 
evidence for a significant extinction (~ 25 per cent) of the A com- 
ponent at 5 GHz is very strong. The relative shift of the source be- 
tween these two frequencies seems to be negligible. This work will 
be extended in the future by comparing medium resolution data sets 
at different frequencies in order to study the propagation effects 
in the lensing galaxy in more detail. Observations with MERLIN 
at different frequencies will allow us to measure the position de- 
pendent Faraday rotation and depolarisation, providing invaluable 
information about the ISM of the lensing galaxy. 

A direct measurement of the galaxy position in B02 18+357 
by HST/ACS observations will be available soon. The formal un- 
certainty of this measurement may be larger than the LensClean 
model constraints but it will be an absolutely independent and 
therefore complementary direct measurement. The comparison will 
either confirm the lens models or give information about deviations 
of the true mass distribution from the relatively simple models. We 
will then reach the point were the uncertainty in Ho is no longer 
dominated by the lens models but by the time-delay uncertainties 
themselves. A reanalysi s of existing monitoring data l Big gs et"all 
1999; Cohen et al. 2000) and new monitoring campaigns will there- 
fore allow further improvements. With all the constraints avail- 
able now or soon and the simplicity of the possible mass models, 
B02 18+357 has the potential to lead to the most robust measure- 
ment of Ho of all time-delay lenses. 
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APPENDIX A: USING THE SUBCOMPONENTS AS 
MODEL CONSTRAINTS 

Applying the lens equation Q for two nearby positions gives to 
first order: 

Az s = M _1 Az (Al) 
Here M denotes the magnification matrix 

(M _1 )ij = Sij - . (A2) 

Clearly any information about the small scale structure of the im- 
ages A and B (flux density ratio, subcomponent positions and 
shapes) can only provide constraints for the relative magnification 
matrix M^ 1 M B . 



Al Relative Positions 

It would be possible to use the subcomponents of A and B as im- 
ages of two independent sources to constrain the models in the 
usual way. Since the separation of the subcomponents is much 
smaller than the typical scale of the system, however, these data 
can only constrain the first derivative of the deflection angle, i.e. the 
magnification matrix. A different approach to constrain this matrix 
directly is therefore more appropriate. 

Given a model value for the relative position vector of r B in 
the source plane, the corresponding vectors in the lens plane are 
given by Eq. tAli . leading to a contribution to \ 2 °f 

X 2 = ( r » -Mir^Cr 1 (n -M 4 f B ) (A3) 

i=A,B 
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with the measured lens plane vectors of ta, J"b and correspond- 
ing covariance matrices Ca and Cb- This equation defines an or- 
dinary linear least squares problem leading to a weighted mean of 
the backprojected ri for the solution: 




MiC^Mi ^MiC^Vi (A4) 



A2 Shapes and sizes 

With the common polar ellipticity parameters, fits can only be done 
numerically and are potentially very unstable. The following Carte- 
sian formalism allows a direct analytical fit for the elliptical shape 
of the source to minimize the deviations from the observed images. 
An ellipse with major and minor axis a and b with direction of the 
major axis (mathematical sense) can be described by the matrix 
E. For all points x located on the ellipse, the following equation 
holds: 

x ] E~ x x = 1 (A5) 

(R+ + R- cos 20 R- sin 20 \ 
~ \ R- sin 20 R+ - R- cos 2(f>) { ' 

R + = ^ (A7) 

2 i2 

R- = ^ (A8) 

The astronomical position angle (north through east) is related by 
p.a. = 7t/2 — <f). Transformation between source and lens plane is 
linear for the ellipticity parameters Eij : 

E s = MX'EaMX 1 (A9) 
= Mb'EbMb 1 (A10) 

We can therefore use the same simple general linear least squares 
formalism as in Appendix I A 1 1 with transformation matrices for the 
three-dimensional ellipticity vector (R+,R- cos 20, R- sin 20). 
In this way the best source plane ellipticity matrix/vector and the 
residual contributions can be calculated analytically. The Carte- 
sian approach also avoids singularities and degeneracies for certain 
shapes, e.g. for circular sources. 

This paper has been typeset from a TjX/ ETjX file prepared by the 
author. 
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